Ising model on a hyperbolic plane with a boundary 
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Abstract 

A hyperbolic plane can be modeled by a structure called the enhanced binary tree. We study 
the ferromagnetic Ising model on top of the enhanced binary tree using the renormalization-group 
analysis in combination with transfer- matrix calculations. We find a reasonable agreement with 
Monte Carlo calculations on the transition point, and the resulting critical exponents suggest the 
mean-field surface critical behavior. 
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Critical phenomena on a plane has been one of the most well-studied area in statistical 
physics. One obvious reason is that it is a plane that we can most easily visualize, but 
it is also because a two-dimensional (2D) system is often analytically tractable. It is true 
that one could say the same about one-dimensional (ID) systems as well, but collective 
orders in a low- dimensional system tend to be fragile against thermal disorders; it is from 
the dimensionality d = 2 for many model systems to exhibit a phase transition at a nonzero 
temperature T > 0. These remarks are best illustrated by the Ising model defined by the 
following Hamiltonian 



where J > is the ferromagnetic interaction strength, and h is the magnetic-field strength. 
The first summation is over every pair of nearest neighbors, and spin at site k can take 
its value from ±1. As is well known, the ID Ising model does not have any magnetic 
order at finite temperatures, while the 2D counterpart undergoes a continuous order-disorder 
transition at a finite coupling strength K = (iJ . Here the inverse temperature is denoted 
as /3 = (ksT)'^, where kg is the Boltzmann constant In this regard, the lower critical 
dimension of the Ising model is 2. 

By planes, however, we do not always have to mean the flat geometry. There are many 
kinds of curved planes observable in physical or biological structures, and critical phenomena 
on such planes may exhibit intriguing features. One can assign the Gaussian curvature to 
a plane, which is positive (negative) when the plane looks like a part of a sphere (saddle). 
With a constant positive Gaussian curvature, the plane will be eventually closed to form 
a sphere, whose radius is inversely proportional to the size of the curvature. This implies 
that the curvature should vanish at every local point if we are to work with a very large 
system size. Even though there can remain global topological constraints determined by 
the positive curvature, many statistical-physical properties will converge to those of the 
fiat geometry in the large-size limit. In the hyperbolic geometry with a negative constant 
curvature, on the other hand, the magnitude of the curvature is not necessarily coupled to 
the system size, and therefore this case is often regarded as more suitable to study effects 
of the curvature. The price is that the surface area A expands exponentially as its radial 
length scale grows. It immediately leads the boundary of the surface dA to expand at the 
same rate, so we find that the boundary fraction dA/A never vanishes even in the large-size 
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limit. The thermodynamic hmit is not uniquely defined for this reason. For example, one 
may use the periodic boundary condition ^ and the behavior will not necessarily be the 
same as with the open boundary condition. By neglecting the presence of the boundary, it 
has been argued by many authors that the Ising model will undergo a mean-field-like phase 
transition on a negatively curved plane: a Monte Carlo analysis away from the boundary 
suggested convergence to the mean- field exponents in Ref. |3i], which was later supported by 
the corner transfer-matrix renormalization-group (RG) method ^. The Ginzburg-Landau 
theory provides a qualitative explanation for this mean-field behavior when only the bulk 
part is considered 5]. For rigorous results under transitivity, one may refer to Ref. Q and 
references therein. 

By accepting this nonvanishing boundary as a part of physics, it becomes possible to study 
nontrivial statistical-physical properties due to the boundary, as has been done in Refs. 5|, 
ffl. A co,„„„,Jy used approach to simulate such a plane is to begiu with hyperbohe 
tessellation with regular polygons and truncate the lattice generation at a certain layer |3I.I10|. 
An advantage is that it makes every point equivalent except at the boundary layer. But 
there is also an alternative lattice structure for such a plane, called the enhanced binary 
tree (EBT) as shown in Fig. [H^a) ll| . It is obtained by adding links between branches in a 
binary tree [Fig. U^b)], so the EBT itself is not a tree, strictly speaking. Although it is not 
a uniform tiling of a hyperbolic plane but made up of triangles and tetragons, this structure 
is actually easier to study analytically: for the bond-percolation problem, for example, it is 
possible to argue that the eniergence of a single giant cluster occurs at a critical occupation 



probability p = Pc = 1/2 12,113!]. This analytic tractability is particularly important because 
numerical calculations suffer from the exponential growth of the system size and therefore 
can only give very rough estimates. There can be in fact one more transition point where 
the correlation diverges in percolation or the Ising model js, 9|, but in this work, we focus 
on the order- disorder transition where the order parameter becomes nonzero. The simple 
binary tree in Fig. IH^b), for example, cannot have the latter type of order- disorder transition 
at any finite temperatures j^. 

This Brief Report is intended to extend the analytic approaches for percolation to the 
Ising model on the EBT. We mainly rely on real-space RG methods and compare the results 
with the numerical data obtained by Monte Carlo (MC) calculations. Combined with the 
transfer-matrix method, this approximate RG calculation allows us to estimate the critical 
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FIG. 1: Schematic descriptions of (a) the EBT and (b) the simple binary tree, drawn on the 
Poincare disks. The number of layers is L = 10 in both the cases, (c) The actual structure used in 
our MC calculations. Here L = 4 is shown as an example, and the structure is extended further 
down for larger L. In (c), we remove the two upper most horizontal connections since they either 
make a self-link (the top dotted line) or double links (the dashed lines) under the periodic boundary 
condition in the horizontal direction. 
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FIG. 2: (Color online) (a) Magnetic order parameter and (b) Binder's cumulant of the Ising model 
on the EBT structures with different sizes. The symbol oo in panel (a) means extrapolated values 
according to Eq. 



point and critical indices as well. The results suggest the mean-field critical behavior of free 
surfaces. 

Let us begin with numerical methods used in this work, and give a very rough estimate of 
the critical coupling K^. of the EBT. For MC calculations, an EBT structure is constructed 
with a certain number of layers, L, as shown in Fig.[T](c), where a periodic boundary condition 
is imposed in the horizontal direction. The system size is then given as = 2^^^ — 1, 
which is an exponential function of L as mentioned above. The Ising Hamiltonian [Eq. ([1])] 
is simulated by the Wolff single-cluster algorithm The magnetic order parameter is 
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FIG. 3: (Color online) (a) Recursion scheme based on the finite-lattice method. The double lines 
represent coarse-grained effective couplings, (b) Iteration by Eq. ([5]) at various K's. The straight 
line has slope 1, i.e., Kn+i = Kn- 



defined by 
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where the bracket (. . .) means the thermal average [Fig. [2](a)]. A convenient quantity to 
locate the critical point is Binder's cumulant defined by 



(|m|^) 



3(|m|2)2' 

whose crossing point suggests Tc ~ 2.45(2) in units of J/fc_B [Fig. EJj^b)]. We extrapolate the 
magnetic order parameter under the assumption that [sl 



(|m|) ~ aN-'^ + b 



(2) 



where the parameters a, b, and are found by the least-square fitting at each T. That is, 
we choose the fitting parameter (p that best describes {\m\) as a linear function of A^"*^. We 
then observe that the limiting value of {\Tn\) at — t- oo could vanish at T > 2.4 [Fig. EJ^a)], 
which supports the above estimation of T^. The exponent is found to have (p ~ 0.2 around 
T 

-I- c- 

In order to study the problem analytically, let us consider the block-spin transforma- 
tion [Ij. In Fig. El^a), we illustrate the transformation in terms of K, which is indexed by 
the iteration step n. The bare coupling K will be hence identified with Kq. This block-spin 
transformation replaces the upper spin block having two triangles and intermediate bonds 
by a single bond, mapping Kn to Kn+i- We consider such a transformation at the outmost 
boundary layer so that an EBT with L layers can be mapped to another EBT with L — 1 



layers. The block-spin transformation can be performed by the majority rule at each triangle 
without ambiguity. The results are 

^g+K„+i ^ g5i<:+3K„ _^ ^^-K+mr, ^-3K+3K„ _^ 2e^^+^" 

_^^-5K+K„ + 2e^-^" + 4e-^-^" + 2e-^-3^" (4) 

with a certain analytic function g, which appears as a consequence of removing short length 
scales. Therefore, one obtains a recursion relation 



1 



Gl{K,Kn) 



(5) 



from Eqs. ([3]) and (jlj). We are interested in a fixed point Kn = Kn+i = K^o- Note that this 
fixed point is always stable, if it exists, since the slope at the crossing is less than 1 [Fig. |3]^b)]. 
So it cannot be a point of phase separation, and we should instead ask ourselves when the 
renormalized coupling strength K^o becomes infinite. It makes sense since it is equivalent to 
asking when the renormalized occupation probability becomes one in the bond-percolation 
problem [l^. The limiting K^o diverges to infinity when K ^ Kc ^ 0.465 810 [Fig. [Ht^b)]. 
However, this estimate corresponds to Tc < 2.2, which is not supported by the MC data 
above. We have also tried working with a larger spin block, but it hardly improves the 
estimate. In spite of such defects, which will be commented on again later, this block-spin 
transformation may give a glimpse of how to perform an RG analysis on this structure. 

In order to find a recursion scheme explaining the MC results better, let us consider using 
transfer matrices. The idea is to describe exactly a single layer first and then replace it by 
a ID chain. It is straightforward to write the transfer matrix for the ID Ising chain with 
coupling strength Kn+i as 

(f,K„+i p—K„+i \ 
K K • W 

g— K„+i gA„+i I 

with two eigenvalues A^""^^^ = e^"+i + e~-'^"+^ = 2 cosh Kn+\ and Aa""^^^ = 6-^"+^ — e~-'^"+i = 
2 sinh i^'n+i On the other hand, the upper spin block in Fig. IH^a) has three spins and 
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FIG. 4: (Color online) (a) Recursion scheme based on the transfer-matrix method. The double 
lines represent coarse-grained effective couplings, (b) Iteration by Eq. ([7|) at various K's. The 
straight line has slope 1, i.e., Kn+i = Kn- 

therefore yields an 8 x 8 transfer matrix T' in terms of K and Kn- The largest eigenvalue is 
obtained as 



1 

— ( 
2 



2e 



with 



RiiK,Kn) = [(e'^"+2)(e'^ + e2^")e=^™" + e^™"+e^^"] 



The second largest eigenvalue is also available as 

2e4i^+2if„ _ ^^2K ^ ^)2g4X„ ^ ^^4K ^ ^yK+6K„ ^ ^R^{K,Kn) 



with 



R2{K,Kn) = e^^+^^"{[-3coshir-cosh3/s: + 3cosh(ir-2ir„) + cosh(3/s: + 2ir„) 
+4 cosh Kn cosh(K + Kn) sinh 2K]'^ - 32 sinh 2K sinh^ 2ir„} . 

Since the matrix T' has six more eigenvalues Ag""* > A4"''' > ■■• > Ag"\ we may need to 



consider eight eigenvalues, or equivalently, eight energy levels from 



(n) 



IuAS") to 



E. 



(n) 



In Ag"'' in total. However, assuming that the first excitation from the ground 



state determines the most dominant behavior, we approximate this with a two-level system 
described by Eq. ([6]). According to this approximation, the first energy gap £'2"^ ~ -^1"'' 



between the two lowest levels should be set equal to i?2"^^'' ~ -^1 



n{n+l) 



In 



7 



in the simphfied two-level picture. In short, we keep the ratio between the two largest 
eigenvalues at each iteration by 

'^2 ^2 



aS") Ai 



(n+l) ' 



which can be compared to Eq. ^ in our first attempt (see Ref. 15|] where a similar idea 
is applied to the quantum Ising chain). It can be also interpreted as adjusting correlation 
lengths since the correlation length ^ from a transfer-matrix calculation is given as 

-1 



(n) 



In 
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Equation ([7]) defines a recursion relation for getting Kn+i out of K and Kn- We depict some 
cases of different K values in Fig. Hl^b). The limiting coupling strength K^o becomes infinite 
when K ^ Kc ^ 0.416 550 7, which corresponds to ~ 2.400 668. Therefore, compared 
to the previous block-spin transformation, this approach describes the MC data better. 

From Eq. (|8]), we can see how the correlation length ^ along the boundary behaves as K 
approaches K^. The behavior turns out to be 

-1 



\K-KA~'/\ 



(9) 



I indicates the mean-field 
and the MC analysis in Refs. js], 16| . 



\n[tanh K^{K)] 

It is related to correlation observed in a sufficiently inner part of the system, but mediated by 
the boundary. The correlation-length exponent u = 1^/2 in Eq 
result, which is consistent with the prediction in Ref. 
We again note that a number of boundary layers are mapped to a ID chain that we are 
looking at, and that Eq. is obtained in this respect. Although traveling along a single 
boundary layer is not the shortest path between a pair of spins, the renormalized boundary 
does contain the shortest path so that the mean-field exponent is found in this RG sense: it 
has been argued that the actual bare correlation will be an exponentially decaying function 
due to the radius of curvature js]. By introducing magnetic field hi at this renormalized 
boundary part, we find 



Ai = 6-^°° cosh (3hi + \l e^^^ sinh^ (3hi + e^^^^ 



which leads to the local susceptibility at /ii = as 

1 d^Xi 



/3Ai dhl 



I3e 



hi=0 



From Koo{K) shown in Eq. ([9]), it is straightforward to find 

Xi--\K-K,\-'/^. (10) 



In the surface critical phenomena jl7|, the correlation decays as G\\{r) ~ r^'^+^^'^ii in the 
parallel direction to the free surface, while it decays as G±{r) ~ ^-d+2-'q^ perpen- 
dicular direction. According to the RG theory, two-spin bulk correlation G{r) ~ ^-d+2-ri 
is related to a scaling dimension by G{r) ~ f""^^*^, where Xh = d — Uh with the RG 
eigenvalue Uh. Suppose that the boundary scaling operator has another scaling dimension 
Xh,s- Then G\\{r) ~ r~'^^'^-'' and G±{r) ~ ^-^h-^^h.s depending on where the spins lie, and it 
is therefore predicted that ri\_ = [r] + ri\\)/2. The local susceptibility diverges as \K — Kd'"^^ 
with 7i = z/(2 — ?7_l). At the upper critical dimension, r^y = 2 and f] = 0, and the exponent 71 
is thus predicted to be 1/2 in the mean- field theory. Therefore, Eq. flTUl) is fully consistent 
with the mean- field description as well as Eq. This mean- field result also gives a clue 
to the poor performance of the block-spin transformation above: the magnetization distri- 
bution gets more and more spread instead of having sharp double peaks around ±1 as the 
dimension increases, so the concept of a block spin becomes less adequate. 

In summary, we have studied the Ising model on top of the EBT structure by using 
the RG methods combined with the transfer-matrix calculation of its single layer. It is 
an approximate calculation considering only the two largest eigenmodes at each iteration, 
and it is hard to assess the error in this approximation, as is usual in many real-space RG 
methods. We have nevertheless reasonably predicted the order- disorder transition point, 
which suggests that this approximate calculation can capture some essential features of 
the system. Our main finding is that the critical behavior predicted from this analysis is 
consistent with the mean-field theory of free surfaces. Since many authors have expected 
that the bulk part far away from the boundary will exhibit the bulk mean-field transition, 
our RG analysis supports the mean-field picture all the way up to the boundary part in 
understanding the Ising model on a hyperbolic plane. 
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